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ABSTRACT 

The quasar ULAS J1120+0641 at redshift z — 7.085 has a highly ionised near zone 
which is smaller than those around quasars of similar luminosity at z ~ 6. The spec- 
trum also exhibits evidence for a damping wing extending redward of the systemic 
Lya redshift. We use radiative transfer simulations in a cosmological context to inves- 
tigate the implications for the ionisation state of the inhomogeneous IGM surrounding 
this quasar. Our simulations show that the transmission profile is consistent with an 
IGM in the vicinity of the quasar with a volume averaged Hi fraction of (/hi)v ^0.1 
and that ULAS J1120+0641 has been bright for 10 6 -10 7 yr. The observed spectrum is 
also consistent with smaller IGM neutral fractions, (Jhi)v ~ 10~ 3 -10~ 4 , if a damped 
Lya system in an otherwise highly ionised IGM lies within 5 proper Mpc of the quasar. 
This is, however, predicted to occur in only ~ 5 per cent of our simulated sight-lines 
for a bright phase of 10 6 -10 7 yr. Unless ULAS J1120+0641 grows during a previous 
optically obscured phase, the low age inferred for the quasar adds to the theoretical 
challenge of forming a 2 x 10 9 M Q black hole at this high redshift. 
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1 INTRODUCTION 

The spectra of quasars with redshifts z > 6 are an important 
probe of the ionisation state of the intergalactic medium 
(IGM) in the early Universe. The intergalactic Lya opac- 
ity rapidly rises with increasing redshift, and consequently 
there has been much debate whether these quasars have 
been observed before the Universe was fully reion ised (e.g. 
iBecker et ai]|2007l : lLidz et al.ll2007l ; lMesingedl201oh . Studies 
of the small windows of Lya transmission observed in close 
proximit y to these objects pro v ide important clues in this 



Loebl 120041; iMesinger fc Haimanl |2004 



respect JWvithe 

Fan et al.ll2006l ; iBolton fc Haehnerdl2007al ; I Alvarez fc Abel 



20071 ; iMaselli et al.ll2009l ). These h ighly ionised "near zones" 
lie between the red edge of the iGunn fc Peterso"nl (|l965l ) 
trough and the quasar Lya emission line, and arise because 
the hydrogen close to these sources is highly ionised. 

iMortlock et al.l (|201ir i (hereafter Mil) have recently 
reported the highest redshift quasar found to date, 
ULAS Jl 120+0641 at z = 7.085, discovered using data 
from the UKIRT I nfrared Deep Sky Survey (UKIDSS, 
lLawrence et al.ll2007h . In their analysis of the Lya near zone 
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observed in a low resolution (R ~ 1400) VLT/FORS2 spec- 
trum of moderate signal-to-noise, Mil find the extent and 
shape of the near zone is distinct from those around the 
much debated z ~ 6 quasars. Assuming the quasar is em- 
bedded in a homogeneous IGM, Mil concluded the observed 
transmission profile of ULAS J1120+0641 is consistent with 
a surrounding IGM for which the volume averaged neutral 
fraction is (/hi)v > 0.1 at z ~ 7.1. 

In this Letter we examine the implications of the Mil 
observation in more detail using simulations of radiative 
transfer through an inhomogeneous IGM. These simulations 
have successfully explained the properties of Lya near zones 
observed at z ~ 6 in the context of a highly ionised IGM 
jBolton et alj|20ld ). We begin in Section 2 by giving a brief 
overview of Lya near zone sizes. In Section 3 we describe 
our numerical simulations before comparing them directly 
to the Mil measurement in Section 4. Finally in Section 5 
we discuss the implications of our results and conclude. We 
assume the cosmological parameters Sl m = 0.26, Qa = 0.74, 
Qbh 2 = 0.023 and h = 0.72 throughout, and unless other- 
wise stated all distances are given in proper units. 
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Figure 1. Summary of near zone size measurements in the cur- 
rent literature. Left: Near zone sizes against the quasar absolute 
magnitude. Bight: Near zone sizes against red shift. The sizes h ave 
been rescaled by a factor io 4 ( 27 + M i4so)/3 jFan et al.ll2006l) to 
correspond to a common absolute magnitude of M1450 = —27. 
Note this scaling assumes -Rnz °c N 1 ^ 3 , where TV is the emission 
rate of ionising photons from the quasar, appropriate if the near- 
zone traces the boundary of the quasar H II region. If the near- 
zone extent is instead set by resonant absorption, -Rnz °c N 1 / 2 
llBolton fc Haehneltjr2007al) . 



2 QUASAR NEAR ZONE SIZES 

At the luminosities typical of observed z ~ 6 quasars, and 
assuming an optically bright phase of t q ~ 10 7 yr, the size 
of a quasar Hn region expanding into a partially neutral 
IGM and that of a proximity zone emb edded in an already 
highly ionised IGM are comparable ijBolton fc Haehneltl 
l2007al ). Further, if the surrounding IGM is significantly neu- 
tral, the red Gunn-Peterson (GP) damping wing will re- 
duce the observed size of the near zone and Lya absorp- 
tion may extend redward of t he quasar systemic redshif t 
ijMiralda-Escude fc Reel Il998l : iMesinger fc Haimanl liooi ). 
Unfortunately, the possibility of confusing this damping 
wing with a collapsed region with a high H 1 column density 
in an otherwise highly ionised IGM complicates matters fur- 
ther. These issues make the interpretation of the Lya near 
zones observed around these q uasars with respect to the 
IGM ionisation state ambiguous llBolton fe Hachnclt 2007al : 
iMaselli et alJl2007l : lLidz et al.ll2007f) .~ 

A summary of the existing measurements of Lya near 
zo ne sizes arou nd high redshift quasars, compiled recently 
bv lCarilli et"aD (|201oh . is displayed in Fig. [1] Interestingly, 
ULAS J1120+0641 has a smalo near zone size given its 
bright absolute AB magnitude, M1450 = —26.6. This is 
demonstrated by the evolution of observed near zone sizes 
with redshift, rescaled to correspond to a common magni- 
tude of M1450 = —27, in the right-hand panel of Fig. Q] 
Mil have noted that the small near zone size around ULAS 
J 1120+0641, coupled with Lya absorption redward of the 
systemic redshift from a putative red damping wing, may 
indicate the IGM surrounding ULAS J1120+0641 is signif- 
icantly more neutral than around quasars at z ~ 6. We 
now investigate this in more detail using realistic synthetic 



The only quasar with a similarly small ne ar zone is CFHQS 
J0210-0456 at z = 6.44 jWillott et al.ll2010h , but this object is 
considerably fainter and is consistent with the other z ~ 6 quasars 
on accounting for its fainter luminosity. 
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Figure 2. Top: Four simple models for the redshift evolution of 
the background Hi photo-ionisation rate along a quasar sight- 
line. The grey shading represents the range covered by the 
Lya forest to rest frame Ly/3 fo r a quasar at z = 7.085 . The 
measurements were obtained by IWvit he &: Boltonl i201lT ) from 
the Lya forest opacity and by ICalverlev et al. I feOllli " using the 
proximity effect. Note that for UVB4 we set Thi = at 2 > 6.9. 
Bottom: The neutral hydrogen fraction averaged over 100 simu- 
lated sight-lines produced by the UVB models including (left) and 
excluding the quasar (right). The averaging is performed for dis- 
play purposes only, and results in neutral regions being displayed 
as spikes with (/hi) ~ 10 -2 . The crosses mark the position of 
all DLAs in the sight-lines constructed using UVB4, and the grey 
curve in the right panel shows the logarithm of (/hi) for UVB1 
excluding the self-shielding approximation, multiplied by a factor 
of 0.01 for clarity. 



absorption spectra which employ a variety of evolutionary 
histories for the IGM ionisation state. 

3 NUMERICAL SIMULATIONS 

We construct Lya absorption spectra using a high resolution 
hydrodynamical simu lation combined wi th a line-of-sight ra- 
diative transfer code (jBolton et al.l 2010). The hydrodynam- 
ical simulation was performed using the GADGET-3 code 
jSpringelll20"05h and has a box size of lOh comoving Mpc 
and a gas particle mass of 9.2 x 10 4 /i~ 1 M o . A total of 100 
sight-lines of length 96.6/i _1 comoving Mpc were extracted 
around haloes identified in the simulation to construct skew- 
ers through the IGM density, temperature and peculiar ve- 
locity field. The largest halo mass in the simulation volume 
at z — 7.1 is 4.7 x 10 10 h~ 1 M ( 7 1 , which likely underestimates 
the ULAS J1120+0641 host halo mass. However, resolving 
the low density IGM is more important for modelling the 
Lya absorption, and a small, high resolution sim ulation pro- 
vides the best compromise jBolton et alj|2010h . Following 
Mil and iBolton et a l. (2010), we assume the spectrum of 
the unobscured, radio quiet quasar ULAS J1120+0641 is a 



broken power law where /„ oc v~ -° 5 for 1050 A < A < 1450 A 
and /„ oc v~ " 15 for A < 1050 A, yielding an ionising photon 
emission rate of N = 1.3 x 10 57 s _1 . 
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Table 1. Average Hi fractions and the incidence of optically thick absorbers within 0.05 Mpc < R < 5Mpc of the quasar host halo in 
all 100 simulated sight-lines including the self-shielding prescription. The absorber columns are estimated by integrating the H I densities 
over a scale of 20 kpc. From left to right, the columns list the volume and mass averaged Hi fraction, the number of LLSs (10 17,2 cm -2 < 
Njji < 10 20 ' 3 cm -2 ) and DLAs (Nm > 10 20 ' 3 cm -2 ) per unit redshift and the fraction of sight-lines where there is at least one LLS or 
DLA, all prior to the quasar turning on. In the final four columns the latter are also given after a quasar with age 10 6 yr and 10 7 yr 
has ionised the surrounding IGM. Note that for models UVB3 and UVB4, where the IGM has a large neutral fraction, identifying LLSs 
becomes ambiguous. At z = 7.085, Az = 1 corresponds to ~ 43.8 Mpc. 
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UVB1 


5.0 x 10" 4 


2.7 x 10~ 2 
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0.3 


0.33 


0.03 


0.16 


0.03 
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0.03 


UVB2 


6.2 x 10~ 3 


8.3 x 10" 2 


55.4 
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0.99 
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0.07 
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UVB3 


1.1 x lO" 1 


3.4 x lO" 1 




0.8 
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0.07 




0.06 


UVB4 


1.0 


1.0 




1.0 




0.11 




0.08 




0.07 



The ionisation state of the IGM along the sight-lines 
prior to the quasar turning on was initialised using four dif- 
ferent models for the background Hi photo-ionisation rate, 
displayed in Fig. [2] However, it is important to note that 
reionisation is an inhomogeneous process, and as a result 
each model may be realised at the same time i n different 
regions of the IGM jMesinger fc Furlanettcl l2008). Further- 
more, the ionisation state of the IGM in the vicinity of a mas- 
sive host halo will be biased due the clustering of low er lumi- 
nosity sources (|Lidz et alj|2007l ; IWvithe et al.ll200sh . These 
models are thus designed to explore a range of different ini- 
tial H i fractions in the IGM close to the quasar, but are not 
representative of global reionisation histories. 

As we only compute the radiative transfer for the quasar 
ionising radiation, we also add a simple prescription for re- 
gions which are self-shielded from the ionising background. 
Assu ming that the t ypical size of an H i absorber is the Jeans 
scale (ISchavdIixHh . the hydrogen number density where 
an Hi absorber begins to self-shield can be approximated 
as n H ~ 3.6 x 10" 3 cm" 3 (Fhi/IO" 12 B - l ) a/a (T/10* K) 2 / 15 . 
We use this expression to compute the density threshold 
at which hydrogen is self-shielded prior to being ionised by 
the quasar, and set the background photo-ionisation rate to 
zero in these regions. The weak temperature dependence is 
computed self-consistently using the gas temperature in our 
simulations. The lower panel of Fig. [5] displays the result- 
ing neutral hydrogen fraction averaged over all 100 sight- 
lines assuming t q — 10 7 yr and including (left panel) and 
excluding (right panel) ionising radiation from the quasar. 
The grey curve in the right panel also shows the H i fraction 
for UVB1 excluding the self-shielding model. Self-shielded 
regions become rapidly more common for ultra-violet back- 
ground (UVB) models with decreasing Hi photo-ionisation 
rates, increasing the volume averaged neutral fraction and 
producing a patchy ionisation structure along each sight- 
line. A summary of the average Hi fractions and incidence 
of optically thick absorbers in all models which include the 
self-shielding prescription is given in Table 1. Note that the 
quasar photo-ionises optically thick absorbers in its vicinity, 
significantly reducing the incidence of these systems. How- 
ever, the impact of this process is reduced for shorter quasar 
ages or the smaller ionising photon mean free path associ- 
ated with a larger IGM neutral fraction. 

4 RESULTS 

We compare the measured transmission profile around 
ULAS J1120+0641, described in detail by Mil, to two of our 



simulated sight-lines in Fig. [3] For each sight-line we con- 
sider quasar ages of 10 6 yr and 10 7 yr, and for each quasar 
age the effect of optically thick absorbers is illustrated by 
models which include and exclude the self-shielding approxi- 
mation. This qualitative comparison highlights an important 
point: for an initially neutral IGM (UVB4, orange curves) 
or a partially neutral IGM (UVB3, blue curves) and a short 
quasar bright phase (10 6 yr) a prominent damping wing with 
absorption extending redward of rest frame Lya is present. 
However, optically thick absorbers will also have an impact 
on the near zone transmission profile redward of Lya if a 
collapsed neutral region is close to the quasar; the presence 
of a damped Lya absorber (DLA) with Nm — 10 20,5 cm -2 
at 3.1 Mpc, marked by the arrow in the right main panel, 
demonstrates this. 

We may consider the relative likelihood of these possi- 
bilities by analysing the transmission profiles for all 100 of 
our simulated spectra. Considering the moderate signal-to- 
noise and low resolution of the observed spectrum of ULAS 
J1120+0641 and the difficulty of continuum fitting due to 
the complexity of the intrinsic Lya and N v emission, we use 
two simple properties to characterise the transmission profile 
for a quantitative comparison. These are: the distance from 
the quasar whe re the transmit ted fraction first falls below 10 
per cent, 7?nz l|Fan et al 1I2006T ), and the transmitted fraction 
at the rest frame Lya wavelength, Ti2i6- These are measured 
after smoothing the spectra with a box car window of width 
2.4 A (rest frame) to remove fluctuations on smaller scales 
due to noise and absorption lines. The corresponding val- 
uefl for ULAS J1120+0641 are R NZ = 1.9 ± 0.1 Mpc and 
T1216 = 0.80 ± 0.07 after removing the intrinsic absorption 
line at rest frame Lya. 

Fig. [4] displays box- whisker and scatter plots of the dis- 
tribution of these two quantities. Both quantities not only 
depend on the assumed UVB model but also on how long the 
quasar has been emitting ionising photons. For example, if 
self-shielding is included and the quasar has been bright for 
only 10 6 yr, both quantities for UVB3 are strongly affected 
by neutral gas in Lyman limit systems (LLSs). Note also 
the extended tail in the distribution of 7?nz and T1216 when 
self-shielding is included for UVB1 and UVB2. This arises 
due to the stochastic occurrence of LLSs and DLAs close to 
the quasar, and is more prominent the shorter the time the 



2 The measurements include a systematic error due to the uncer- 
tain continuum placement (see Figure 4 in Mil). 
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Figure 3. Simulated quasar spectra constructed using the four different UVB models displayed in Fig. [2] The spectra are processed to 
mimic the Mil spectrum by convolving with a Gaussian instrument profile with FWHM= 240 kms -1 , rebinning onto pixels of width 0.4A 
(rest frame) and adding Gaussian distributed noise using the observed noise array. The observed transmission of ULAS J1120+0641 is 
shown by the light grey curve in each panel. Left: Sight-line 1. Note the impact of the red GP damping wing on the spectrum for models 
with a fully neutral IGM (orange curves). Right: Sight-line 2. In this sight-line a prominent red damping wing is seen for all UVB models 
when the self-shielding prescription is included. A DLA with TVjji — 10 20 5 cm -2 is present at the position marked by the arrow. 



quasar has been bright and thus able to photo-ionise proxi- 
mate optically thick absorbers. Fig. [4] suggests that the near 
zone around ULAS J1120+0641 is due to either (i) a par- 
tially (fully) neutral surrounding IGM with (/hi)v ~ 0.1 
(1.0) and a quasar bright phase of ~ 10 6 yr (10 7 yr) or (ii) 
the presence of a high column density Hi system in close 
proximity to the quasar. It is important to note, however, 
that (i) and (ii) are not mutually exclusive. In our modelling 
an abundance of optically thick absorbers close to a quasar 
is a generic prediction due to the shorter ionising photon 
mean free path expected when the IGM is only partially 
ionised. Our simulations confirm the suggestion of Mil that 
the spectrum of ULAS J1120+0641 is consistent with an al- 
ready very highly ionised IGM ({/hi)v ~ 10" 3 -10" 4 ) if a 
DLA lies close to the quasar. However, this occurs in only 
~ 5 per cent of our sight-lines for t q = 10 6 -10 7 yr. A fully 
neutral IGM for t q — 10 6 yr is ruled out by the near zone 
size, but is consistent with the observed transmission for 
a bright phase duration of 10 7 yr. However, ionisation by 
lower luminosity sourc es expected to c l uster around a mas - 
sive quasar host halo ijLidz et al.ll2007l : IWvithe et al.ll2008l ) 
may make this extreme possibility less likely in practice. 

We should caution, however, that if we have underes- 
timated the number of proximate DLAs present at z ~ 7, 
if the ionising spectrum of the quasar is much softer than 
assumed here, or if the environment of the quasar host halo 
is s ignificantly overdense out to 2 Mpc from the quasar (but 
see lCalverlev et al ] 20 111 ) .then we will overestimate the like- 
lihood of a significantly neutral surrounding IGM and/or 
underestimate the age of the quasar. 

5 DISCUSSION AND CONCLUSIONS 

We have presented a large suite of realistic synthetic trans- 
mission spectra which model a variety of evolutionary his- 



tories for the ionisation state of the IGM surrounding the 
quasar ULAS J1120+0641. The transmission profiles depend 
on the age of the quasar, the photo-ionisation rate in the 
surrounding IGM and the presence of high column density 
absorbers arising from collapsed, high density regions. For 
the same initial conditions for the IGM ionisation state the 
Lya transmission profile varies greatly due to the stochas- 
ticity of the optically thick absorbers and IGM density dis- 
tribution along different sight-lines. 

Inferences regarding the ionisation state of the 
surrounding IGM for Lya transmission around ULAS 
J1120+0641 are correlated with the duration of the optically 
bright phase of the quasar. For a short bright phase dura- 
tion of 10 6 yr, the data are consistent either with an IGM 
with a volume averaged neutral fraction (/hi)v ~ 0.1 or a 
highly ionised IGM with (/hi)v ~ 10" 3 -10" 4 and a proxi- 
mate DLA. On the other hand, a bright phase of 10 7 yr is 
consistent with both a highly ionised IGM and proximate 
DLA or a fully neutral IGM. However, we find only ~ 5 per 
cent of our simulated sight-lines exhibit a proximate DLA 
capable of producing a red damping wing, implying the for- 
mer scenario is less likely (but see also the caveats in section 
4). In either case, this suggests the quasar has not been 
bright for significantly more than ~ 10 7 yr; a longer bright 
phase will produce a near-zone which is too large, either by 
enlarging the quasar Hn region or photo-ionising progres- 
sively more proximate optically thick absorbers. A timescale 
of 10 7 yr is still a factor of five less than the e-folding time 
of a black hole growing at the Eddington rate and radiating 
with 10 per cent efficiency, as might be expected for a moder- 
ately spinning black hole. For a rapidly spinning black hol e 
the discrepancy would be even larger (jSiiacki et al.ll2009l ). 
Unless the black hole grows in an earlier optically obscured 
phase, this adds to the challenge of building up a black hole 
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Figure 4. Left: Box-whisker plots of the Lya near zone size (upper panels) and transmitted fraction at the Lya rest frame wavelength 
(lower panels) measured from synthetic quasar spectra. For each UVB model the distributions arc shown with (SS) and without the self- 
shielding prescription. The thick horizontal line in each box displays the median of the distribution, and the vertical extent of the broad 
and narrow boxes correspond to 68 and 95 per cent of the distribution around the median, respectively. The whiskers display the range 
of the simulated data. The horizontal dashed line and grey shaded region in each panel displays the measurement for ULAS J1120+0641. 
Right: Scatter plots of the transmission at the Lya rest frame wavelength, T1216, against near zone size, i?NZ- The data point with error 
bars corresponds to ULAS J1120+0641. 



with an estimated mass of 2 x 10 9 M Q so early in the history 
of the Universe. A higher quality spectrum will hopefully ex- 
clude or confirm the possible presence of a proximate DLA 
by enabling the identification of associated metal lines and 
aid in narrowing constraints on the neutral fraction of the 
surrounding IGM. 

More general conclusions regarding the ionisation state 
of the IGM at z ~ 7 are speculative due to the large fluc- 
tuations in the ionisation state expected on scales much 
larger than probed by the ne ar zone of ULAS J1120+0641 
(Mcsi nger fc Furlane tto 2008). If the neutral fraction in the 
general IGM at z ~ 7 is indeed (/hi)v ~ 0.1 this is neverthe- 
less consistent with a constant or weakly decreasing ionising 
emissivity between z ~ 6 — 7. The rapid change in the photo- 
ionisation rate by two orders of magnitude or more in our 
models is balanced by a similar change in the number of 
LLSs and thus the m ean free path of ionising photons {e.g. 
iMcQuinn et~ai]|201ll ). Since a constant or perhaps even ris- 
ing emiss ivity is required for the completion of reionisatio n 
at 2>6 (|Miralda-Escudell200o1 ; rBolton fc Haehneltj|2007rj ). 
this may be compatible with the pos sibility that rei onisation 
was yet to fully complete by z — 6 (Mcsinger 2010]). Lastly, 
the possibility of a significant neutral fraction in the IGM 
surrounding ULAS J1120+0641 make it an excellent target 
for searching for extended Lya emission from the expanding 
ionisation front (|Cantalupo et al1l2008l ) as well as studies of 
the IGM ionisation state using the redshifted 21 cm transi- 
tion. 
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